Gravity from Discrete Interactions, 
a laboratory of ideas for field theory 



Manoelito M de Souza 
- - - Universidade Estadual do Sudoeste da Bahia * 

^ . Departamento de Ciencias Exatas - Vitoria da Conquista - BA - Brazil 

O '■ and 
^ I Universidade Federal do Espirito Santo* 

^ ■ Departamento de Fisica - Vitoria - ES - Brazil 

Q: 
^: 

Abstract 

O ■ The observed accelerated expansion of the universe is an indication that somehow, in some 

conditions, gravity may become a repulsive interaction. The very concept of discrete interac- 
^ I tions, compatible with General Relativity as an effective theory of averaged fields, can handle 

this and other pressing problems without the need of searching for ad hoc exotic sources. Semi- 
quantum models of discrete interactions based on quantum exchanges between point-like masses 
are a rich laboratory for enlightening field interactions and classical-to-quantum transitions in 
^ ■ field theory. A very simple version is presented for exhibiting how naturally it makes contact 

' with inflation, dark matter and dark energy issues. 
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1 Introduction 

In few more years the scientific community will be celebrating the centennial year of the two great 
^SJ , pillars of modern theoretical physics. Quantum Mechanics and General Relativity, without being 

able to make them compatible, harmonize them into a single scheme. Quantum Mechanics carries 
the idea of interaction discreteness realized through exchanges of a quantum interaction, a packet 
of energy and momentum, while General Relativity associates the gravitational interaction to the 
smooth geometry of a curved spacetime. This incompatibility is a serious problem because they 
■ are both unbeatable on fitting the experimental data with high accuracy on their respective domain 

of validity and applicability. General Relativity overseeds the old Newtonian scheme reproducing 
its dynamics in its Newtonian limit, a non-relativistic weak field limit. There are great challenges 
facing General Relativity [T]. The observed accelerated expansion of the universe is an indication 
that somehow, in some conditions, gravity may become a repulsive interaction. The old question 
of the flat rotation curves of stars and galaxies is still unsettled. These questions have in common 
very large scales of distances which lead to the above mentioned weak field Newtonian limit In 
the Newtonian scheme there were no room for a repulsive gravity but General Relativity, assuming 
its validity on such scales of distance, can handle it in various ways but none without controversy. 
The Einstein field equations are the equality between, in one side, the geometric one representing 
the observational data, the field, the test mass acceleration, and, in the other side, the sources, the 
matter content represented by the stress-energy tensor. Feeding this scheme with the observed data 
for finding the matter content, an exotic fluid that fulfill this data, may be an easy but certainly 
dangerous procedure as it raises the Popper's falsifiability question if there is no pre-established 
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^The flat rotation curves occurs when the field is smaller than 10~^m/ which, by the way, is of the same order 
of the Pioneer anomaly [2]. This may not be a mere coincidence. 



criteria on the acceptability of what is found. The same observations are vaUd for the galactic 
flat rotation curves. Then we have the dark matter and the dark energy controversies. Inflation 
in cosmology, its nature, what powers it and what smoothly disconnects it is another controversy, 
or unsettled question in gravitation theory. These are certainly among the greatest riddles for the 
modern researcher in gravitation theory. 

The objective of this letter is to point to a new remarkable fact that they can altogether be 
enlightened with just very simple ideas on discrete interactions and no further else assumption. For 
discrete interaction it is meant the instantaneous emission/absorption of a quantum of interaction 
triggered by the instantaneous absorption/emission of a quantum of interaction, with free, straight- 
line propagation between any two consecutive interactions, in a flat spacetime background. As an 
illustration we present here the simplest discrete interaction model for gravity between two point- 
like masses in a non-relativistic relative radial motion; it can easily be associated to a cosmological 
model of a spherically symmetric distribution of dust. Such kind of naive models, with no free 
parameters and great variety of levels of complexity, are useful laboratories of ideas in the search 
for a better understanding of field interactions and quantum to classical transitions. The resulting 
dynamics, strongly dependent on the distance scales and on the initial conditions, oscillates between 
repulsive and attractive phases. It starts with a repulsive exponential expansion that smoothly 
changes into a Newtonian phase to return later, again to a repulsive expansion. This rich scenario 
needs no extra input to develop besides the two point-like masses and discrete interactions. The 
well known success of General Relativity on high accurately fitting some experimental data does not 
stand against gravity being a discrete interaction since the very General Relativity i3 , like Classical 
Electrodynamics, can be seen as an effective theory of averaged interactions in a discrete interaction 
context [4l[5l|6]. They are both finite discrete interaction theories; singularities and other well known 
problems are shown to be consequences of the averaging process in the field definition. 

2 The model 

Let us consider two point-like masses, M and m, in relative non-relativistic radial motion under a 
discrete gravitational interaction, m is a test mass (Af >> m) and its momentum and position 
(relative to M) are given by 

ri-l 

i=o 

as the momentum only changes at discrete points when there is the emission/ absorption of an 
interaction quantum, a graviton, and 

r„ = r„_i + Ar„_i = r„_i + ^^M^-i = + V (2) 

m ^ — ' m 

3=0 

as there is free straight-line propagation between any two consecutive interactions. The label n 
counts the number of discrete interactions starting from at the initial position rg with the ini- 
tial momentum p^. With such discrete interactions a trajectory is piecewise-linear with vertices at 
the interaction events, where velocity, as the trajectory tangent vector, is not defined because the 
discrete interaction causes a sudden change of velocity. So the notion of acceleration, if not as an 
approximation, is meaningless since it is null in the free propagation period between two consecutive 
discrete interaction and is not defined at the vertices. 

The kinematics 

Among many other alternatives, we take for our model of classical discrete interaction 



(3) 



where a;„ = ^ and a is the kinematic coefficient that describes the system free propagation between 
any two consecutive interactions. This choice is justified assuming that, let us say, the graviton, 
propagates with the speed of hght c and, crudely, that the emission of an interaction quantum is 
triggered without delay by the absorption of an interaction quantum. Then, in the non-relativistic 
approximation 

At„«^ = ^x„, (4) 
c c c 

The dynamics 

The dynamics for our model will be chosen with an eye on the above mentioned non-relativistic 
weak field limit, which, after (jS]), implies on 



as the most immediate choice, since 



Apn^P/Xn, (5) 

Ai„ dt 

n 



requires 

/3 1 _ GMm GMm 1 



(7) 



which fixes the dynamic coefficient /3 as 



p^J^. (8) 
cro 

There is no free-to-be-adjusted constant in this toy model of discrete classical interactions between 
two point-like masses. However, we strongly remark that ([5|) is not the unique possible choice, 
even after ([3]), and neither the most appropriate considering that the Newton's equation is just an 
approximation to a correct formulation of the problem. It is, nonetheless, interesting to see that 
discrete interaction can produce repulsive gravity even in such a quasi- Newtonian simple context. 

On the other hand, ([6]) has a price as it requires an approximate continuity on t„ and p„ (grossly 
Ai„ << 1 and Ap„ << 1) which, with ([3]) and ([5]), leads to 

/? « Xn « -, (9) 
a 

that, for dimensional reasons, must be seen just as a qualitative indication of the existence of such 
a range of validity for the Newton's law of gravity. 

A word of caution 

A word of caution, with respect to is however necessary. We are dealing here with the grav- 
itational interaction between two point-like mass. In practice this would imply on something like 
the gravitational interaction between two electrons, or neutrinos or some other elementary particle. 
For a long time still to come this would be a purely academic business! If M is due to a number 
Nm of point-like masses distributed on a region of radius vm, with tm << tq, and uncorrelatedly 
interacting with m, then a = ^ must be replaced by 

- (10) 



cNm 



Just for an order of magnitude, we remind that light takes about 8 minutes to reach Earth from the 
Sun which cannot, obviously, be the order of the Ai„. A^j\/ for the Sun is such a huge number that 
most certainly we will never be able to measure such a tiny lapse of time At„. On the other hand, 
in this non-relativistic approach a is the only place where c, the speed of light, will appear. We will 



keep N]\i = 1 everywhere and in the case we need to consider a real physical observation involving 
a composed object we just have to replace c by cNm- 

c^cNm (11) 

In general we do not know Nm but we can assume, in a gross approximation, that 

Nm oc M. 

A similar reasoning [5^ must be considered at the moment that m absorbs an interaction quantum 
and if m is composed of point-like objects then the dynamic coefficient (3 (or the probability that 
this absorption occurs) is proportional to Nm and then to m. This justifies part of the Newton's 
statement on the gravitational force (being proportional to the product of the masses). 



3 Solving for the equation of movement 

After m) and (El), with a = ^ = and b = = ^ (H may be written as 

n-1 j-1 n i-1 ^ 2 

x„ = 1 + — ^(po + J2 >^ = 1 + II — (12) 

u i=o ji=0 j=0 ji=0 ^1 

which can be solved by re-iterations but it is more convenient to use an expansion in series powers 
of a 

Xn — ^ ^ Xn,s-j Xn,Q — 1? (l'^) 
s>0 

where Xn,s are coefficients to be determined in terms of the initial conditions rg, po and of n. We 
find [5] 



Xn,o = 1, a;„4 = a( ^ ) , Xn.k = for /c > [n/2] 



and for k < [n/2]: 



xr.au = BA^-' ( j^J + 2(fc - l)hBA^-^ (^fc - I j + • ■ • ^^^^ 

xn,2k+i = aA'^ (^2fc + ij + + " (v2fcj ^ ■ ■ ■ ' ^^^^ 

where the dots indicate decreasing sequences of n-binomial numbers. We have introduced, for 
convenience, the two constants A and B, defined as 



2_ / _ GMm\ _vl~ w| _ ^2 



^-«^ + 26^£,h^-^ .:^L-^ = ^, (16) 



B = A-b=^^, 



where vq and vea are, respectively, the velocity and the escape velocity at the initial position rg in 
the standard Newtonian potential; Wq = Uq — w|n ^ . yl is clearly related to the system initial energy. 



4 The n-dominant order interaction. 



If we consider in each Xn,s only the contributions from the term with the highest power in n, 
(n >> 1), 

J^ + Oin'-') (17) 
we have 

,2fe 



n 

(2fc) 

^2fe+l 

(2fcTT)! 

and then, from ([T5)) . assuming A ^ 0, 



Xn,2k-A''-^Bj^^ for k>0 (18) 



Xn,2k+i^aA''-^—— for k>0, (19) 



[n/2] 2fr ["/2: 



„2fe t'V^J „2fc+l 



fc>l ^ fc>0 ^ ^ 



with Xn = Q!\/C4n = ^^n, n » 1, and avC? = << 1. It is a good approximation to replace 
the partial sums by their respective asymptotic series, 

x„ w — + — coshX„ H — ^ sinhX„ (21) 
A A J A 



GM 

— (cosh(^„-l), (22) 



with 0„ = X„ + 00 - ^ + (/-o and 

a^/A ^vqWq 



tanh(/)o = 



Replacing in ([T]) the sum by an integral, which is proved to be valid [9] at this order of approximation, 
leads to 

q=Pn - kn ^ Po - ko, with fc„ = 2mwo(e'^" - 1)^"^. (23) 

q, which has the physical dimension of linear momentum, is a conserved quantity as it does not change 
with n. It would be kind of peculiar to have linear momentum conservation under gravitational 
interaction but actually it boils down 9^ to energy conservation 



q = mw = TOW V, 



2GM / , 2GM 



I = mJvl = V2rnE. (24) 



The dominant-order contributions reproduce the usual Newtonian picture of continuous interactions 
with its known interaction potential. 



The "potential energy" 



The "potential energy" , retrieved from the q conservation 

jj ^ <f -yj ^ (g-p«)(g + p«) ^ fcn(2g + fc„) 

" 2m 2m 2m 

can be seen from a new perspective: we can be interested on knowing the zeroes of the potential, 
2q + fcfl = and kn = 0. In this dominant order contribution approximation each condition equally 
leads, with (|23p . to the same point at infinity 

VE,n = 0, or, equivalently 1/^,-1 = 0, 

as expected once the potential is the old Newtonian one, 

m^{vl - «| „) ~pI _ mi>| „ _ GMm 
Zm 2 r„ 

But it is remarkable the possibility of having two distinct zeroes (this can happen with the inclusion 
of sub-dominant contributions) and then there is a place for a repulsive gravity too, as we will see 
next. 



5 Sub-dominant contributions. 

It is convenient now, for wq 7^ 0, to combine the initial condition parameters wo, va [va > 0) and tq 
(or ve,o) into a single one h, h = 



"0 



> 1, if Vq > VEn 



''-^ 1 >0, ifvo<VE,o 



As /i > 1 is finite with the equality to 1 being only asymptotically reached for 

ve,o —^0 or ro 00, 

we will be assuming from now on ft, >> 1 except when the opposite is explicitly stated. 

Including now the next n-order contributions i.e. the terms with the second highest power of n 
in each x„ 



n\_n^- (^) 



kj (fc) 



+ 0(?i*-^) (26) 



we have, instead of ([Ti)) and ((T5|) . 



and 



= BA^-^ + 2ik - l)bBA^-'j^^ + . . . 

" A y {2ky. (2/c-l)! 2 (2fc-2)! j^'" 
Xn,2k+i = aA^ \ ' + ahA^-^ {kA + {k - 1)5) 



(2fc+ 1)! J y • y 11 {^2k)\ 

a ( (nVI)^^+i bB {n^/Af' nB^ {nVAf''-^ \ 
73 I (2A: + 1)! "If (2fc)! + 2A {2k - 1)\ j ' ' ' 



and they lead, with (IT51) . to 



where 



x„ = To + Ti cosh X„ + T2 sinh X„ + ix„M^(r3 cosh X„ + sinh X„) (27) 

r— VP, bB Wn , 

he A^fA 2c ^ ^ 
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- W = 


2hc 


~ h^) > 0, 


T3 = 



2h^ 

> 



T4 = ^ = -/i:; « ft > 

1 - 



The term ^ = ~lf ^~2t«^ ~ marks the sub dominant contributions; being proportional 

to ^ it shows the importance of the velocity of propagation of the signal (quanta) and suggests that 
the next order of interaction may require a proper relativistic approach. On the other hand ^ << 1 
does not assure |^^| < 1, as long as there is no a priori constraint on the size of h. 

Higher n-order contributions 

The coefficients Tq, Ti and T2 carry contributions from both the dominant and the sub dominant 
n-orders and they will receive contributions from any other lower n-order that may be added. Con- 
sidering the possibility of including the N first n-orders, (j27p can be generalized to 

x„ =ro(iV)+^^^^^|i^(r2fc+i(7V)coshX„ + r2fc+2(^)sinhX„). (28) 

Every new n-order added adds contributions to every T-coefficient, including the Tq one. It is a 
peculiar aspect of these n-order expansions that the relevance of lower n-order terms increases as n 
grows in contraposition to the usual approximation expansions in terms of small coupling constants 
where each order of contribution is fixed. This is unexpected as in general the gravitational effects 
decrease with the distance; these non-Newtonian effects increase. As the distance keeps increasing, 
contributions from lower n-orders become more and more relevant. Probably there is a closed form, 
still to be found, for each T-coefficient as well as for the entire equation of which the expansion ([28|) 
would be just a partial sum. On the other hand there is no point on adding new contributions if 
they only become relevant with distances that go beyond our observation horizon. 



6 Discrete versus continuous 

Let us compare (|2T|) and ([27| . They represent two distinct perspectives, two different approaches to 
the dynamics of gravity. (|2T|) stands for the non-relativistic weak-field limit of General Relativity, 
for the vision of an interaction field as continuous as a smooth surface while (l?7)) shows additional 
consequences of accumulated effects of successive discrete interactions. Let us now point to some 
striking differences between them. 

While ([2T|) describes just one phase evolution, ([27|) describes a multi-phase one. The dynamics is 
very sensitive to the initial conditions. \W\ « 1 reduces (f27l) to (|2T|) as it turns the sub dominant 



n-order contributions negligible. On the other hand as n (or equivalently the distance or the time) 
increases, the last two terms on the right-hand side of (|27l) . proportional to WXn, and that we 
are calling the sub dominant contributions, increases (WXn times)faster and may compete or even 
dominate, radically changing the dynamics. Let 

\W\Xn « 1 (29) 

mark the point where these two contributions grow even; for |W|X„ > 1 the non-Newtonian effects 
increase faster than the Newtonian ones but for |M^|X„ << 1 (l27l) reduces to 

Xn^To + Ti coshX„ + Ta sinhX„ for \W\Xn « 1 (30) 

which is just (I21[) with different constants. The dynamics is strongly dependent on the initial 
conditions and, even after this simplification, it still encompass a great variety of physical situations 
of which we mention here some few. The same procedure applied to (j2ip leads [5] now to a distinct 
conservation law 

P , (e"^" - Z+) f3 , (1 ~ Z+) 

(e-" - Z_ ) = - 1^ (13^ 

with 



Z± (Ti +T2) = -To± + r| - T2 = -To±W (31) 
and leads to a dominant logarithmic contribution to the potential 

For 

« a;„ (33) 
both q and [/„ reduce, in a first order approximation, to the Newtonian results, as expected. 

Gravitation between elementary particles 

The interesting result comes from the other extreme condition 

\W\ » Xn (34) 

that reduces ([32]) to 

(/.»-^l„(l + HE) ,35, 

maWVA \ Xn J 

whose gradient produces an effective "acceleration" 

rodUn qGM 1 



m dxn Wwq r„ 

An r„-dependence like that could explain the flat rotation curves of galaxies. Besides, considering 
the conservative potential we can imagine, just to simplify, that a given star was carried to its actual 
state from an energetically equivalent initial position Vq and an initial velocity vq ~ 0, then the 
condition reduces to 

or 

2c2 

On, AT >> 

^0 



where an,N = GM/r'^ is the Newtonian value for the star acceleration at r„ which must be bigger 
than a reference value in a kind of twisted Milgron's condition [8], in reference to an empirically 
found condition satisfied by most of flat rotation curves of galaxies. The nagging point is that 
plugging back the Ni\i factor of dTTI) turns the condition (IMl) into 

(t7 ) >Xn (36 

Nmc 

which can be realized only with few-components objects like elementary particles as they have small 
values of Nm ( Nm = 1), a;„ « 1, excluding any macroscopic object, not to mention a star or galaxy. 
It is saying that the gravitational interaction between two elementary particles could, in principle (it 
depends on h, the initial conditions as the dynamics is very sensitive to them) radically differ from 
the one we know between macroscopic objects in usual conditions! Had we only a dream technology 
capable of measuring such discrepancies! 



7 Discrete interaction and cosmology 

The discrete interaction toy model considered in this paper can easily be associated to a Newtonian 
cosmological model of a spherically symmetrical distribution of dust of which m is the mass of the 
point dust at r„ and M is the total mass of the dust spherically distributed inside the sphere of 
radius r„. It is mostly a consensus that inflation, an early exponential expansion of the universe, 
is a necessary ingredient of big-bang cosmologies despite the many questions lingering on, waiting 
for consistent answers. Questions like, how did inflation started? powered by what? and how did it 
smoothly ended into the present known state of the universe? Similar questions can also be made 
with respect to the observed accelerated expansion of the universe, that seems to require an evolution 
system with a repulsive gravitation phase. Discrete interaction, we believe, have something to add 
on these issues. 

Exponential expansion phases 

In the dynamics described by (1271) there are two main exponential expansion phases, an early and a 
late one. Let us start re-writing it as 

Xn =To+ (ti +T2 + XnW^'' + T4\e " 



2/2 



n-T2+ X„W^^^^ j ^ (37) 
from which we can see that it describes an exponential expansion as long as 

Ti + T2 + X„T4^^i-tIl » Ti - T2 + X„W^^i^ (38) 



Xn« — 4r (39) 



\W\T, 



4 



and also when 



_ 2iTi+T2) ^, 2Ti _ 2T2 

>> |,,,| ^„ — , „ , ~ njrim 



|M^|(r3-|-T4) \W\Ti \W\Ti 

as Ti >> T2 and T4 » T3. It is natural that if we want to use this discrete interaction model 
for pursuing a more comprehensive understanding of gravity we must consider the whole (|27p and 
not just (jSO]) . To repeat with (|27)) the same treatment given to ([2T|) and to ((30)) is not as easy but 
we can do, for now, with numerical computation for a glance on its dynamics. See the figures below 
where we have used the command 
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Figure 1; 3 phases from an infinite periodical series of repulsion, attraction and repulsion. 



Manipulate[Plot[f{x, h), {x, Xmin, Xmax}], {h, h^,,,, hj^ax}] 

of the software Mathematica 8.0 which is very useful as it allows the observation of the dynamics 
dependence on the initial conditions through the explicit variation of h. For the Figure 1, showing 
the repulsive, attractive and again repulsive sequence of phases, we have plotted the "acceleration" 
Vndvn / dxn against X„ for p7p with A < and with the initial conditions partially fixed by taking 
vq — 1000to/s and leaving ro free, hidden in the parameter h which was allowed to vary in the 
range from 10"'^ to 10^ and x (that stands here for X„) from to 2tt. Increasing the range of x 
would mean including distances beyond, perhaps, our observation horizon; for consistency this would 
require including more n-order contributions. A smaller h implies, for a fixed initial velocity, a larger 
escape velocity which means, for a fixed M, a smaller initial position, while a larger h implies ve,o 
closer to vq- For comparison, using (pij) instead of (P7)) . for the same initial conditions, produces 
Figure 2 and an always negative Newtonian "acceleration" 

Vndvn/dxn = "4.5 X 10^(1 + 1 / h^) (cos0.5X„ + hsinO.bXny^ « 
-4.5 X lO^x-^ < 

in its unique attractive phase. 

We should observe that in both figures the scales of Xn in the horizontal axis are not linear 
with the scale of distance x„; they are related by (I^TI) and 1^7^, respectively. In the notation of 
Mathematica , with the chosen initial condition, we write (|27p as 

x„ = ((1 + h^)/2 + 10-5(1 - h'^)/2) + 

((1 - h'^)/2 - 10-5(1 - h'^)/2 - hlO-^x{l - h'^)/2) cos[x] + 

(h - 10-5(1 - h'^)/2 + xlO-%1 - h^f/A) sm[x]. 

Changing h changes the overall aspect of the plot but not the existence of changes of sign in the 
"acceleration" as it is naturally dependent on the initial conditions. The choice of the snapshot 
taken at h — 360 just reflects the author's esthetic bias without any further deep physical concern. 
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Figure 2: The unique Newtonian phase: a„ < 0. 

8 Conclusions 

The objective of this paper of pointing to the relevance of discrete interaction models as laboratory 
of ideas in the search for a better comprehension of field interaction and of its discrete to continuous 
or quantum to classical transitions has been accomplished. Although being based on a very simple 
classical model of discrete interactions for gravity and being mirrored on the non-relativistic weak 
field limit of General Relativity which is just an approximation to a correct formulation of the 
problem it makes contact with actual major issues in the physics of gravity. Most remarkable, of 
course, is the natural appearing of repulsive gravity in a simple quasi-Newtonian context. The next 
question is how far can it go. 
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